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Abstract: We study the twist-2 double-spin asymmetries for the 2-jets and large-p^ hadron 
(vr, A etc) production in semi- inclusive deep inelastic scattering to 0{as) in perturbative 
QCD. After deriving the complete set of the polarized cross section which is differential with 
respect to the transverse momentum, we discuss characteristic features of the azymuthal 
spin asymmetries, using existing parton densities and fragmentation functions at COMPASS 
and EIC energies. 
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1 Introduction 



High energy experiment with polarized beams and targets has opened a new window for 
revealing QCD dynamics and hadron structures. Ongoing RHIC-SPIN, HERMES and 
COMPASS experiments are going to provide us with a variety of data disclosing spin 
distributions inside the nucleon. The planed Electron Ion Collider (EIC) at BNL is expected 
to be a more powerful and sensitive tool for the QCD spin physics. 

In this paper, we study the hadron production in polarized semi-inclusive deep inelastic 
scattering (SIDIS) off proton. The cross section formula with integrated transverse mo- 
mentum of the final hadron, pr, has been derived in the leading order (LO) [1] and in the 
next-to-leading order (NLO) [2] and has been applied to predict spin asymmetries for po- 
larized A production in SIDIS [3]. We extend these studies to derive cross section formula 
for the complete set of the polarized processes which is differential in and is suited, in 
particular, for the large-p^- hadron production. This study should be useful to get further 
insight into the parton densities and fragmentation functions from HERMES, COMPASS 
and future EIC experiment. It is also complementary to studies on large-p^ hadron pro- 
duction in pp collisions at RHIC and HERA-iV [4, 5, 6]. As typical SIDIS processes, we 
study pion (or 2-jets) and A hyperon production. In the leading twist-2 level, the complete 
set of the processes is 



(i) 




e' + nipr) + X or 


e' + 2 jets + X, 


(ii) 


e+p ^ 


e' + A{pT) + X, 




(iii) 


e+p^ — 


>e' + A\pT)+X, 




(iv) 


e + p 


e' + 7r(pT) + ^ or 


e' + 2 jets + X, 


(v) 


e+p 


e' + A(pT)+X, 





where we have used the notation p and A for longitudinally polarized, and p^ and A^^ for 
transversely polarized hadrons, and we restrict ourselves to electron (or muon) scattering 
off proton. In order that a final hadron carries large transverse momentum, another parton 
has to be emitted in the opposite direction. This is an 0{as) effect in perturbative QCD 
and was investigated in [7, 8] for the unpolarized SIDIS with electron and neutrino beams. 
Here we extend the analysis to the polarized cases shown above. For small-pT- production 
in SIDIS, the origin of p^ may be ascribed to the nonperturbative intrinsic transverse 
momentum of partons confined in or fragementing into a hadron. Derivation of the cross 
section formula based on this idea was performed in [9], which is complementary to the 
present study. 

As will be shown below, our cross section formula for (i)-(v) applicable in the large-p^- 
region diverges as 1/p'^ as pt — ^ 0, which is a manifestation of soft and collinear divergence. 
For the pr-integrated cross section, the soft divergence is canceled by the soft divergence in 
the virtual correction diagrams and the collinear one is factorized into the parton densities. 
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Accordingly, 0{as) cross section at pr ~ is to be interpreted as "distribution" with 
respect to the variable pt- In order to get a finite p^-difFerential cross section, one needs 
to resum the effect of soft gluon radiation as was studied in [10, 11]. Derivation of the 
resummed cross section for a low pt region is beyond the scope of the present study and 
will be presented elesewhere. 

This paper is organized as follows. In section 2 we present the formalism to calculate 
the cross sections for (1.1) following [8]. Section 3 presents the final analytic formula for 
the cross section. In Section 4, we present a numerical estimate of the spin asymmetry 
and discuss its characteristic features, using existing parton distribution and fragmentation 
functions in the literature. Section 5 is devoted to summary and conclusion. 

2 Formalism 

In this paper we are interested in the cross section for the process, 

e{k) + A{pA, Sa) ^ e'{k') + B{pb, Sb) + X. (2.1) 

We define 5 Lorentz invariants for specifying the kinematics. The center of mass energy 
Sep for the initial electron and the proton A is 

Sep = {PA + ~ 2pA ■ K (2.2) 
where we ignored the masses. Conventional DIS variables 

x,j = ^ -q' = -{k - k')\ (2.3) 
^Pa ■ Q 

are determined by observing the final electron. For the description of kinematics of the 
final hadron B, we introduce 

= (2.4) 

Pa ■ q 

and the "transverse" component of q which is orthogonal to pa and ps. 



a _ u PB-q a PA-q ^„ 

it — q Pa 

Pa-Pb Pa- 

qt is a space-like vector and we define its magnitude as 



?r -q'- ^^P\ - ^^P%. (2.5) 

Pa -Pb Pa- Pb 



Qt = ■ (2-6) 
Differential cross section for the semi-inclusive production (2.1) can be written as 
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where we use the covariant normahzation for each state and 1/25'ep is the initial flux. L^i^ 
is the leptonic tensor defined as 



(2.8) 



with our convention for the anti-symmetric tensor 60123 = — e°^^^ = 1. The leading twist-2 
contribution to the hadronic tensor W^'^ can be written as 

fl (ix dz r - — - 
W^^ipA^PB^q) = — ^Tr MA{X,PA,SA)MB{Z,PB,SB)H^''{XPA:PB/Z,q) 



x2Tid [ixpA + q-pB/zr) , 



(2.9) 



where Ma{x,pa, Sa) and Mb{z,pb, Sb) are the distribution and fragmentation functions for 
the hadrons A and B, respectively, and H^^'^{xpa-,Pb/ z.q) is the corresponding hard part. 
Tr indicates the trace over relevant spinor or Lorentz indices. Since we are interested in the 
twist-2 cross sections, pa and pb can be regarded as light-like. To define the complete set 
of distribution and fragmentation functions, we introduce another set of light-like vectors 
n {n? — 0) and w [w"^ — 0) for pa and pb, respectively, by the relation pa • n — 1 and 
Pb ■ w = 1. The complete set of the twist-2 quark distribution for the proton A is defined 
as [12, 13] 

MI,.{x,pa,Sa) = J —e''^{pASA\^jmi{^n)\pASA) 

11 1 
= 2^Aq{x) + 2^aJ5^aM{x) + -;pb$A±i>A^q{x) H — , (2-10) 

where the transverse spin vector S^^ = (0, Sa±) satisfies Sa± ■ Pa = ^ and X\ + S\j_ = 1 
with the longitudinal component of the spin vector, and -I — • stands for twist-3 or 
higher. q{x), A.q{x) and 5q{x) are, respectively, spin-average, longitudinally polarized and 
transversely polarized (transversity) quark distribution. Sq is chiral-odd and hence does 
not mix with gluon distributions. The gluon distrbition is also defined as [12, 14, 15] 

dX 



MT^{x,pa,Sa) = ^ / —e'^-{pASA\tr n,G'^"(0)n,G'-'^(An) \paSa) 



x^ J 27r 

~G{x)gf_, - ^AG{x)te-^''^pA,n, + ■■■, 



(2.11) 



where is the gluon 's field strength, tr means the trace over the color index for C", 
g'^^j_ — g"-^ —p'\n^ —n'^p^j^. G{x) and AG{x) are, respectively, spin-average and longitudinally 
polarized gluon distributions. Notice there is no transversely polarized gluon distribution. 
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Similarly we define twist-2 quark and gluon fragmenation functions for the spin- 1/2 hadron 
B[12, 16]: 



= + -\bi4b^^{z) + -l^$B±i>B^^{z) + • • • , (2.12) 



M^^-^{z,pb,Sb) = / ^e-^'/^OHw,G^-mBX){BX\w^G'^^{Xwm 



= - AG(z)ie^^'"'pBpW, + ■■■, (2.13) 

where the transverse spin vector S'^_|_ = (0, Sb±) in the quark distribution satisfies Sb±-Pb = 
and X%+Sb± = 1, and g'^^j_ = g'^^ —p%w^ —w'^p^B- Physical meaning of each fragmentation 
function is in parallel with the distribution functions defined above, and we specify them 
by putting " ^ " on the corresponding distributions. 

For the actual calculation of the cross section, it is convenient to work in the hadron 
frame [8] . In this frame and take 

g'' = (0,0,0,-g), (2.14) 

Choosing the a;2;-plane as the hadron plane, pb can be written as 

„M _ ^ ^ , i: ^ ^ - 1^ f2 16) 

In order to write the lepton momentum in the hadron frame, we need to introduce the 
azymuthal angle between the hadron plane and the lepton plane in the hadron frame. 
Then the lepton momentum can be parametrized as 

k'^ = (cosh ■0,sinh'0 cos 0,sinh'0 sin 0, —1) (2-17) 



and k'^ ^ k" - qi^ with 



coshV'^^^^-l. (2.18) 

Q 



With these definitions, the cross section for (2.1) can be expressed in terms of Sgp, x^j, Q^, 
Zf, q^ and in the hadron frame. Obviously, 4> is invariant under boost in the q direction, 
so that the in the hadron frame is the same, for example, in the center-of-mass system 



5 



of the virtual photon and the initial proton A. The expression for these variables in terms 
of the lab. frame variables is given in the Appendix. 

Let us proceed to calculate the cross section. As is shown in (2.16) the transverse 
momentum of Pb is pr = ZfQT in the hadron frame, and we shall derive the cross section 
formula differential in q^. Lowest order diagram in the QCD coupling relevant for the 
processes (1.1) is shown in Fig. 1. We follow the method of [8] for the calculation. To 
isolate the 0-dependence in the cross section, we first expand the hadronic tensor W^" in 
terms of the complete set of the independent tensors. To this end we introduce the following 
4 vectors which are orthogonal to each other [8] : 

T^-^{q' + 2x,,p%), 
f 

Q 



= (2.19) 



These vectors take the form of = (1,0,0,0), = (0,1,0,0), ^ = (0,0,1,0), Z^' = 
(0, 0, 0, 1) in the hadron frame, which greatly facilitates the actual calculation presented 
below. Note that T, X and Z are vectors, while Y is an axial vector. Since W'^'^ satisfies 
the current conservation, q^W^'^ = q^W^^ = 0, it is easy to see that W^'^ can be expanded 
in terms of 9 independent tensors, for which we employ the following: 

Vr = AT^X'^ + F^F^ V^" = g^"" + Z^'Z", V^"" = T^X*^ + X^'T", 

Vr = X^'X'' - VY", Vr = iiT^X" - AT^T'^), Vr = iiX^'Y" - y^AT'^), 

V^^ = {{T^'Y" - Y^'T"), vr = Tf'Y" + y^T^ V^'" = X^'Y" + Y^'X''. 

(2.20) 

The expansion coefficients of W^'^ in terms of these tensors is easily obtained by the inverse 
tensors Vfe for Vfe {k — 1, - • ■ 9) as 

W^''^ = E vr , (2.21) 



k=l 



where 



~ I „ ^ 1 

vr = -{2T^T'' + X^X" + y^y), vr = T^T'', Vr = — {T^X" + X^T"), 

vr = I {X^'X" - Y^Y") , vr = ^ {T^X" - X^'T'') , Vr = ^ i^^y - Y/'X") , 

Zi Zi Zi 
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2 2 2 

(2.22) 

Note that Vji'^ {k = 1, ■■ - ,4,8,9) are symmetric between ^ and u, and the rest are anti- 
symmetric. V^'^ {k = 1, • • • , 5) are even under parity, while the rest are odd under parity. 
For the processes with unpolarized electrons, i.e., (i)-(iii) in (1.1), both L^^" and W^^ are 
symmetric tensors, accordingly V\^^ {k — l,---,4) contributes to W^". For (iv) and (v) 
both L^'^ and W^'^ are anti-symmetric pseudo-tensors due to a 75 in the trace, accordingly 
only Vj!" (k = 6,7) contribute to W'^'^. Actual calculation of Wp^jV^^ in (2.21) involves 
trace over many 7-matrices, which can be carried out conveniently in the hadron frame by 
using the Tracer program. 

With these V^" , we can decompose the 0-dependence of the cross sections by the con- 
traction with L'"' in (2.8). Define Ak (A; = 1, • • • , 9) as 



A = ^ (2-23) 



then one obtains 



^1 = 


1 -I- cosh^V) 




-2, 


A = 


— cos0sinh 2'ip, 


Aa = 


cos 20sinh^'0, 


Aq = 


—2 cosh?/). 


A7 = 


2cos0sinh'i/). 



(2.24) 

From (2.24), the cross section for (i),(ii) and (iii) have ^-dependence of 1, cos0 and cos 20, 
while (iv) and (v) have only 1 and cos dependence. 

For the case of (iii) e + ^ e' + + X, we need to parametrize the spin vector Sa± 
and Sb± which lie in the plane orthogonal to pa and pb, respectively. We write 

Sa± = cos ^ aX^ + sin ^aY" = (0, cos sin 0), (2.25) 

S^^ = cos Qb cos ^bX^ + sin - sin 9^ cos ^bZ^" 

= (0, cos©BCOS$B,sin$B, -sin©BCOS$B) (2.26) 

where the second equality in each equation follows only in the hadron frame. In the hadron 
frame, ^a,b are the azymuthal angles of Sa,b± measured from the hadron plane, and ©b is 
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the polar angle oi pb measured from the Z-axis which has the Lorentz invariant expression 
In the Appendix, we present the expression for $^ and $s in terms of lab. frame variables. 

3 Analytic formula for cross section 

By applying the method described in the previous section, we finally obtain the cross section 
for the large pr hadron production in DIS in the following form: 



dxbjdQ'^dzfdq^dcf) SirxljS^pQ'^ ^ Jxmin x Jzf 

where Ak is defined in (2.24), ae — e^/An is the QED couphng constant, and we introduced 
the variables 

x^% z^'-l, (3.2) 
X z 



and 



For a given Sep, and qt, the kinematic constraint for x^j and 2;/ is 



'-'ep 



For a given Q'^, Xhj and Zf, qt can take 

<qT <Q 



(3.6) 



Roughly speaking, < < Q corresponds to the current fragmentation region and qr > Q 
corresponds to the target fragmentation region. As is shown in (2.16), the transverse 
momentum of in the hadron frame is 



PT = ZfQT < ZfQ, 
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In (3.1), [/ -D (g) CTjfc] should read the following form for each process in (1.1) (Cp — 
4/3): 

(i) e+p^e' + {7r,A} + X[7]: 

+ j:elq{x)G{z)ar, (3.8) 



where 



at" = 2Cfxz 



1 / Q' 



(g^-g|)')+6} 



af = 2af = 8Cfx?, 

a|^ = 4C^x2-i-(Q^ + g^), (3.9) 

Q(lT 

af = x(l-x)|^(4:^-4 + 2) + 10-^-n, 

af = 2af = 8x(i-x), 

= 2af = 8Cfx(l-?), 

sr.4c.s(i-.f^{e'^+^}. (3.11) 

(ii) e+p^e' + A + X: 

[/ ® D ® afc] = E Ag(a:) Ag(;2) A^af + ^ e;AG(a:) Ag(;^) A^af 

+ 5:e^Ag(x)AG(^)A^af, (3.12) 



where 

A^ar = af (fc= 1,2,3,4), 

^ (2£ - 1) {Q\x - If - q^x^} 
Q^qj'x{x — 1) 

A^af = A^af = 0, 

2{Q^(x-l)-g^x} 

^^^^ - (^^V^^ g- 

(in) e+p^ ^ e' + A^ + X: 

[/ ® ^fe] = ^ el5q{x)5q{z)/:^TaV. 

where 

Ayaf = 4CfC0s($a-$b), AT^f^O, 

Qt qr 

(iv) e+p^ e' + {n,A} + X: 

[f0D®a,] = Y: elAq{x)q{z)/\LodV + E AG(a;)g(z) A^o?f 

+ EejAg(a;)G(^)A^o?f, 
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where 



Q(lT 
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2x-l( x-lQ^\ 
. _2g(^-l)(2£-l)^ 

(v) e + p ^ e' + A + X: 

[f®D®ak]=J2 elq{x)Aq{z)AoLar + E ejG(x) Ag(z) AoL^f 

+ Ee5g(a;)AG(;.)AoLar, (3.22) 

where 

AoLae'j = A^oaej, (3-23) 

2x2-2x + l/ Q2X 
Ao.a- = ^ 

. 2Qi^-mx-l)^ 
Qt z 

Several comments are in order here: 

(1) From the above formula one can separate the 0-dependence of the g^-differential cross 
section for (i)-(v) in (1.1) as 

with (72 = for (iv) and (v). Transverse momentum of partons also gives the same (j) 
dependence as a purely kinematic effect [17]. It is, however, interesting to test whether 
the above 0(0^5) formula gives quantitatively right magnitude for each component. 
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(2) The cross section depends on Sep through Ak in (2.24). For XbjSep/Q'^ ^ 1 (see 
(2.18)), we have ^1,3,4 ^ Aqj (0-dependence factored out), which results in the 
strong suppression of the asymmetry for (iv) and (v) compared to (ii) and (iii) in 
(1.1). 

(3) In the large-gr SIDIS, contribution from the gluon distribution and fragmentation 
functions is of the same 0{as) effect as the quark contibution, so that g-r-differential 
cross section is expected to be a more sentitive tool to determine detailed form of the 
polarized gluon contibution than the gr-integrated case [8] . 

(4) Except for (iii), (Tq has a nonintegrable l/g|i- dependence at qx — ^ 0, while cri,2 are 
integrable. For the case of transverse polarization (iii), a2 has the l/g|,-dependence, 
while (To,i is integrable at qx —>■ 0. To get a finite meaningful result for all value of 
qx, one needs to include higher order effect based on the resummation technique as 
in the case of the spin-average case [10, 11]. 



4 Numerical estimate 

4.1 Azymuthal asymmetry 

In this section, we present a numerical estimate of the polarized cross sections (ii)-(v) in 
(1.1) in comparison with the spin averaged one (i), using the exisiting parton distributions 
for the nucleon and the fragmentation function for pion and A. For this purpose we consider 
the azymuthal asymmetries defined as follows: 

, > _ Jo dQ'^dxbjdzfdq^dcj) _ (Jq 

\ ISaSb = r2n 11 dSg-'"' ~ ^av ' 

JO "■y^dQ2dxijdzfdq^d4> 

I / iw JU V V ^ / dQ^dx},^dzfdqt,d(p (7 n\ / ^ -i \ 

(cos(n</.))5,5, = .2.,^ = {n = 1, 2) (4.1) 

Jo "-V dQ2dxtjdzfdq^d<l> 



where a"''" and 0"^°' are the spin-averaged and polarized cross sections and their decompo- 
sition o-Q^i^2^ is defined in (3.26). The subscript SaSb in the left-hand-side of (4.1) specify 
the spin states of the hadrons A and B for each process in (1.1). We use the symbol O 
for spin-average, L for longitudinal polarization and T for transverse polarization. Accord- 
ingly, SaSb in (4.1) can be SaSb — OO, LL, TT, LO, OL corresponding to (i), (ii), 
(iii), (iv) and (v), respectively. By definition (l)oo = 1- For all cases, |(l)5^5s| < 1 and 
|(cos(n0))5^s^| < 1. 

The asymmetries (4.1) are still functions of 5 variables, Sep, Q^, Xhj, Zf and qx- We 
estimate them at COMPASS {Sep = 300 GeV^) and EIC {Sep = 10^ GeV^) energies with 
typical kinematic variables where our perturbative formula are valid. Comparison of pre- 
dicted asymmetries with experimental data may be better achieved by integrating over some 
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of the variables to gain statistics. Here we will not try this procedure, but will show the 
nonintegrated azymuthal asymmetry, intending to show a typical dependence of the asym- 
metry on each variable. To determine the kinematic variables, we first note from (2.18) 
and (2.24) that at large coshi/j = 2xbjSep/Q'^ — 1, the spin asymmetry for (iv) and (v) is 
strongly suppressed compared with (ii) and (iii), so that a smaller value of cosh^^ is prefered 
to get large asymmetry. In addition, if one chooses the same values of and cosh^ at 
COMPASS and ETC energies, variation of the asymmetry at different Sep is wholly ascribed 
to the variation of distribution and/or fragmentation functions and the partonic hard cross 
sections at different x^j. Keeping this in mind we have chosen (Q'^jX^j) = (100 GeV^,0.4) 
for the COMPASS energy and {Q^Xbj) = (100 GeV^ 0.012) for the EIC energy, by which 
"valence" and "sea" region of the parton densities are mainly probed. 

As a reference distribution and fragmentation functions, we use GRV parton density 
for the unpolarizcd nucleon[18], GRSV polarized parton density [19], KKP fragmentation 
function for vr^+vr^ [20] and FSV A-fragmentation function [3] for polarized and unpolarized 
A + A. In all cases we use NLO versions, putting the factorization scale /i^p = Q^. 

4.2 Pion or 2-jets production 

The process (iv) e+p e' + 2 jets + X or e' + 7r + X may be useful to get further insight into 
the polarized parton distribution in the nucleon. Here we present an estimate of (1)lo and 
{cos 4>)lo for the charged pion production. Figure 2 shows the gr-dependence ((a) and (b)) 
and ^/-dependence ((c) and (d)) of these asymmetries at COMPASS ((a) and (c)) and EIC 
((b) and (d)) kinematics, {cos (f))oo is also shown for comparison. For the spin asymmetry, 
we showed the results with standard (solid line) and valence (dashed line) scenarios in the 
GRSV distribution. At Xbj — 0.4, the spin asymmetry (1)lo can be as large as 50 %, while 
(cos0)/,o is within a few % level. At x^j = 0.012, these magnitudes are somewhat reduced, 
since the polarized parton density is more suppressed in the small x-region compared to the 
unpolarized parton density. At x^j — 0.4, the two scenarios give totally different behavior 
for (1)lOj which can be, in principle, a sensitive test for the polarized parton density. 
Typical magnitude of {cos 4>)oo is larger than {cos 4>)lo as expected, since the unpolarized 
parton density is larger than the polarized one. 

4.3 A production 

Figure 3 shows the ^T-dependence ((a) and (b)) and the ^/-dependence ((c) and (d)) of the 
azymuthal spin asymmetry for (v) e+p — >• e' + A + X ((a) and (c) for {1)ol, and (b) and (d) 
for (cos0)ol) at the COMPASS kinematics. Figure 4 shows the same asymmetries as Fig. 
3 but for the EIC kinematics. In [3], three sets of polarized fragmentation function for A-|-A 
have been constructed based on three different assumptions at a low energy input scale: 
Scenario 1 corresponds to the nonrelativstic quark model picture, where only polarized s- 
quark is assumed to fragment into the polarized A, i.e., AS = Ad = 0. Scenario 2 is based 
on the assumption that the A fragmentation function has a similar flavor decomposition as 
the A distribution function constructed from proton's structure function gi by the SU{3) 
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flavor symmetry, i.e., —Au — —Ad — 0.2 As. Scenario 3 assumes three flavors contribute 
equally to the fragmentation process, i.e., Au — Ad — As. In all three cases. As has 
positive sign. We calculated the asymmetries with these three sets. In both kinematics, 
(l)oL can be as large as 50 %, while (cos 0)o£ is within a few % level, varying with different 
scenarios for the A fragmentation functions. 

The pattern of the calculated asymmetry (1)ol can be easily understood [3]. At the 
COMPASS kinematics with Xf,j — 0.4, contribution to the asymmetry is basically from the 
valence w-quark distribution (see (3.1) and (3.3)). Thus scenario 1 in which only s-quark 
fragments into A gives negligible asymmetry, while scenarios 2 and 3 gives, respectively, 
negative and positive asymmetries, the former being smaller in magnitude. At the EIC 
kinematics with smaller xi,j {x^j = 0.012), a large contribution comes from the sea-quark 
distribution, which is nearly flavor symmetric. Accordingly scenario 1 also gives a positive 
asymmetry with its magnitude smaller than the scenario 3, since only s-quark fragments into 
A. In scenario 2, u and d quark contribution is partially canceled by s-quark contribution, 
resulting in the small negative asymmetry. 

Figures 5 and 6 show the asymmetry for (ii) e + p^e' + A + Xa.t the COMPASS 
and EIC kinematics, respectively. We show the results with the standard scenario for 
the GRSV polarized parton density. For comparison unpolarized azymuthal asymmetry 
(cos (j))oo and (cos20)oo are also shown in Fig. 5, the former being essentially the same as 
the one for the pion production shown in Fig. 2 as expected. Typical magnitude of 
(cos0)ll and (cos 20) at x^j = 0.4 is a few 10 %, a few %, and O(10~^). At Xbj = 0.4, 
and (cos^)^^ have approximately the same magnitude as (1)ol and (cos(0))ol- At 
Xf,j — 0.012, (1)ll and {cos(f)) ll are, respectively, smaller than (1)ol and {cos (f))oL, since 
polarized parton density is more suppressed at smaller x region compared with unpolarized 
parton density. The distinction among three scenarios for the polarized A fragmentation 
function can be understood in the same way as the asymmetry for (v) e + p^e' + A + X 
shown in Figs. 3 and 4. 

Figures 7 and 8 show the spin asymmetry for (iii) e+p^ — > e'-|- A^ +X at the COMPASS 
and EIC kinematics. The calculation is done at cos($^ — $b) = I. In the present estimate, 
we assume the transversity distribution and the trasversity fragmentation function are equal 
to the longitudinally polarized distribution and fragmentation function, respectively: We 
put 5q = Aq and 5q = Aq. This assumption may be justified at a low energy scale. 
But even in that case, they should be different at high energy, since their Q^-evolution is 
different, in particular, at small a; [21]. One can alternatively estimate the upper bound 
based on the Soffer's equality as was done for + p A^ + X [5]. Here we intend to 
present a simplest estimate to study characteristic features of the asymmetry by the above 
ansatz. One sees from the figures that (1)tt is of a few % ~ 20 % at xi,j — 0.4 depending 
on the scenarios for the A fragmentation function. This value is smaller than (1)ll, which 
is expected since there is no gluon contribution to the transverse spin asymmetry. On the 
other hand, (cos0)tt is from a few to 10 %, and even (cos20)tt is of a few % level, which 
are much larger than (cosn(p)LL {n = 1,2). The same feature persists also at x^j = 0.012, 
although absolute magnitude of the asymmetry becomes smaller by the same reason for 
(cosn0)ii {n — 0, 1, 2). It is interesting to check these peculiar features of the hard cross 
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section in future experiments and to get information about transversity. 

5 Summary and conclusion 

In this paper, wc have studied the double spin asymmeties for 2-jets and large-p^ hadron 
production in semi-inclusive DIS with one-photon exchange. The cross section formula for 
the complete set of the spin dependent processes (1.1) has been derived to 0{as). Separating 
the dependence on the azymuthal angle (p between the hadron and lepton planes, the formula 
can be written as 

^ (7o + cos 0-1 + cos 2(f) (72, (5.1) 

where the last term is absent for the process with polarized lepton beams. This decompo- 
sition is purely kinematic, while the magnitude of each component depends on the origin 
of the transverse momentum of the final hadron. At low pt, nonperturbative component 
such as intrinsic-A;^ of partons is expected to contribute as studied in [9]. At moderate 
Pt, higher order effect represented by the resummed cross section is presumably important. 
Our 0(0:5) formula should be useful to clarify the importance of these effects and to see 
at which pt the experimentally observed spin asymmetry follows the perturbative QCD 
prediction. 

In order to see the qualitative behavior of each component in (5.1), we calculated the 
azymuthal asymmetries using the parton distribution for the nucleon and the fragmentation 
functions for vr and A at x^j = 0.4 and Xiij = 0.012. We found that the asymmetry is very 
sensitive to different scenarios of parton distribution and fragmentation functions as was 
expected from the previous studies on SIDIS and pp collisions [3, 4, 5]. Typical magnitude 
of asymmetry from each component in (5.1) turns out to be 0(0.1 ~ 0.6), O(10~^) and 
O(10~^) from the first to third term in (5.1). For the process (iii) e + p^ ^ e' + A'^ + X in 
(1.1), however, the second and third contribution in (5.1) is significantly larger, while the 
first term is smaller, compared with other spin asymmetries. 

As in the case of unpolarized cross sections, ctq component of the polarized cross sections 
have l/g|.-singularity at low qx, except for the above process (iii) where (72 has the same 
singularity. To get a meaningful formula at low qt, we need resummation technique [10, 11], 
which will be reported in a future publication. 
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Appendix 

In this appendix, we give the expression for the Lorentz invariants and hadron frame vari- 
ables in terms of the variables in the Lab. frame. (See [8] for the detail.) In the Lab. frame. 
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4- vectors which appeared in the text can be parametrized as 

p'X = Ea{i,o,o,i), 

= £^B(l,sin^BCOS0B,sin^Bsin0B,cos^B), 
F = 0,0,-1), 
k'^ = £;'(l,-sin^,0,-cos^), 

= {E - E' , E' sme,Q, E' cose - E), 
>S^^ = (0,cos$i,sin$i,0), 
^^j^ = (0, cos6'b COS0S cos — sin sin 

cos sin 05 cos + cos0ssin$£, — sin cos $^), (A.l) 

where $^ ^ is the azymuthal angle of the transverse spin vector Sa,b± which satisfy Sa± ■ 

Pa = Sb± ■ Pb — 0- In order to give the expression for q^, we introduce the polar angle 9^ 
for in the Born amplitude (diagram obtained by removing the gluon line from Fig. 1(a)), 
which allows us to write pb oc (1, sin 0, cos 9^). It is easy to show 9^ is given as 

cot ^ = (A.2) 

With this 9^, we have 

2 -4^'(2^-E')(l + COS^) r . 2... . . . . . 2/i /^n1/aoN 

qr = ^ {sin' {{Ob - 9,)/2) + sin 9b sin sin' (0b/2) } (A.3) 

i — cos u B 

and 



cos = - — ( 1 — 



2gT V XbjS, 



Finally, the azymuthal angles cos$a and cos$b in the hadron frame can be written as 

cos$A = -Sa^ - X ^ — f — £'Bsin6'scos(0s - <l>^) - sin 6' cos ) , (A.5) 

qr \Zf J 



cos $B = ■ „ Sb± ■ Z 
sm v3b 



q2 _|_ q2 

"^Q^ — ~ {e' sin 9 (cos (j)B cos cos — sin 0b sin 

-{E'cos9-E)sm9BCOs^^}. (A.6) 
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! Jc!K 



Figure 1: 0{as) cut diagrams contributing to large pr hadron production in SIDIS. Distri- 
bution function is located in the lower part of each diagram, and <S> show the fragmentation 
function insertion. 
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Figure 2: Azymuthal spin asymmetry for the pion production e-fTr — > e'-fTT-fX. (a) and (b) 
show the ^T-dependence of the asymmetry at COMPASS and EIC energies, respectively. 

(c) and (d) show the 2;j-dependence of the asymmetry at COMPASS and EIC energies, 
respectively. Solid- and dashed- lines, respectively, denote those obtained with standard- 
and valence- scenarios of GRSV polarized parton density Dash-dot lines denotes Azymuthal 
asymmetry (cos 0)oo for the unpolarized cross section. 
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Figure 3: Azymuthal spin asymmetry for e + p ^ e' + A + X at the COMPASS energy, (a) 
and (b) show the gT-dependence of (1)ol and (cos 0)ol, respectively, (c) and (d) show the 
^/-dependence of (1)ol and (cos0)ol, respectively. 
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Figure 5: Azymuthal spin asymmetry for e + p — e' + A + X at the COMPASS energy. 
(a),(b) and (c) show the g-r- dependence of (cos^)^^ and (cos20)ij^, respectively. 

(d),(e) and (f) show the ^/-dependence of (cos0)ll and {cos2(j))iL, respectively. 

Unpolarized azymuthal asymmetry (cos 0)oo is shown in (b) and (e), and (cos20)oo is 
shown in (c) and (f) for comparison. 
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Figure 7: Azymuthal spin asymmetry for e + ^ e' + A^^ + X at the COMPASS energy. 
(a),(b) and (c) show the gfr-dependence of (l)rT, (cos0)Tr and (cos20)rT, respectively. 
(d),(e) and (f) show the ^/-dependence of (1)tt, {cos (f))TT and {cos2(f))TT, respectively. 
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Figure 8: The same as Fig. 7 but for EIC energy. 
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